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Abstract 

We revisit a scenario in which the cosmological constant is cancelled by the potential energy 



' of a slowly evolving scalar field, or "cosmon". The cosmon's evolution is tied to the cosmological 

m ' 

■ constant by a feedback mechanism. This feedback is achieved by an unconventional coupling of 



the cosmon field to the Ricci curvature scalar. The solutions show that the effective cosmologi- 



o 
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cal constant evolves approximately as t~ 2 and remains always of the same order as the density 
' °f ordinary matter and radiation. Newton's constant varies on cosmological time scales, with 

a: 

Gn /Gn *C 1/t. Gn could have been somewhat different, and possibly smaller, at the time of Big 

<D ! 

r~ | , Bang nucleosynthesis. 
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1. INTRODUCTION 

Many ideas have been suggested for solving the cosmological constant problem [lj, but 
none so far has proven to be completely satisfactory An old idea is that the vacuum energy 
of ordinary fields is cancelled by some compensating field, sometimes called a "cosmon" . The 
cosmon would roll down a potential hill until it had lowered the total vacuum energy to an 
acceptable level 0, EH^f- The question is how the cosmon field would "know" when to stop 
rolling. If the cosmon potential happened to have a minimum at exactly the right height, 
it can cancel off other contributions to A. However, that would obviously be nothing other 
than a manual setting of the total A to zero at the minimum, i.e. just the old fine-tuning of 
A that one is trying to avoid. 

If a conspiracy of parameters is to be avoided, it is imperative that the cosmon action 
does not directly depend on the contributions to A coming from other sectors of the theory. 
But how, in that case, is the cosmon to know when to stop rolling? Some sort of feedback 



mechanism is required. Several such mechanisms have been proposed over the years 



y, |5| . The one we will explore in this paper was proposed in the papers of Ref . [2j, y] in the 
mid 1980s, and is based on the idea that the cosmon field couples to the scalar curvature 
R. Since the trace of the Einstein field equations tells us that —R = 2 A + ~(p — 3p), the 
cosmon can learn about A through its effect on R. The conventional wisdom, based on a 
"no-go theorem" stated in Weinberg's well-known 1989 review of the cosmological constant 
problem, is that this kind of feedback mechanism cannot work. However, this "theorem", 
which originally comes from the papers , does not apply to the kind of model we discuss 
in this paper, as we shall see. (Note: we will be using units where WttGn = 1, signature 
-+++, and the sign convention in which a sphere has negative R.) 

Before exploring this idea further, there is a general point that may be important. In 
any realistic feedback mechanism, it would seem that the "effective vacuum energy" A — by 
which we mean the vacuum energy from other sectors of the theory (A) plus the potential 
energy in the cosmon field (V cosmon ) — must fall to zero as time passes in such a way that it 
is always of roughly the same order as the matter energy density, p ma tter- By "matter" here 
we mean all forms of particles, including baryons, photons, neutrinos, dark matter particles, 
etc. The reason for this is as follows. If A falls off significantly more slowly than Pmatter, then 
it will quickly come to dominate the universe, which will end up with f2 maiier being negligibly 
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small. If, on the other hand, A falls off significantly faster than p ma tter, then very early in the 
history of the universe the effective cosmological constant would be much smaller than T 4 , 
where T is the temperature of the universe. However, that would mean that any subsequent 
phase transition (such as one associated with the symmetry breaking of grand unification, 
the weak interactions, or the chiral symmetry of the quarks) would drive A very negative, 
since it would lower A by an amount of order T 4 . Given the arguments of Coleman and De 
Luccia j?J and more recently from a string theoretic viewpoint, Banks 0, it is questionable 
whether such a transition could occur. But if it did, once the universe has negative A, it is 
not clear how it could get out of that hole, as it would have to in order to be consistent with 
present observations. 

The upshot is that any successful feedback mechanism would seem to require that the 
hidden energy density of A plus "cosmon" field energy must be roughly of order the density 
of ordinary matter at all times. This is very suggestive, given the puzzling fact that a dark 
energy has been observed that is comparable to the matter energy today. 

The feedback model we discuss does indeed have the desired feature that A ~ Pmatter at 
all times. As we shall see, however, this does not necessarily account for the dark energy 
that has been observed. The problem is that any mechanism that is designed to erase A will 
tend to erase dark energy too. In the model we discuss there is indeed energy that is "dark" 
(namely A = A + V cosmon ), which generally falls to zero roughly as t~ 2 and has w = 1/3. 
However, after a phase transition A tends to fall more slowly for a period of time, so that w 
can be smaller than 1/3 and perhaps even close to — 1. It is possible that this fact can be 
used to account for the nearly constant dark energy that has recently been observed. 

To return to the specific idea that the feedback occurs through coupling of the cosmon to 
the scalar curvature R, the obvious way for this to happen is through the cosmon potential. 
To illustrate how this might work, consider a toy model with V cosmon = -a(pR } where 

always denotes the cosmon in this paper. The effective vacuum energy, as we defined 
it above, is then A = A — a<pR. As the cosmon rolls down its linear potential hill, A will 
decrease; and as that happens, the scalar curvature will decrease as well. This decrease of the 
scalar curvature reduces the slope (—otR) of the cosmon potential, and therefore the cosmon 
rolling will decelerate (since there is friction due to the Hubble expansion). If we assume 
that R receives contributions only from A, then as A approaches zero, the slope of the cosmon 
potential will approach zero also, and the velocity of the cosmon's rolling will approach zero. 
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In other words, whatever the value of A is, the quantity A will asymptotically approach zero 
as time passes. That is, the potential energy of the cosmon, —a<pR will progressively erase 
A. 

In fact, this is what actually would happen in this toy model, as an explicit solution of 
the field equations shows. Unfortunately, however, Newton's constant suffers a disaster. 
Instead of the usual Einstein-Hilbert action —R, one has —(1 — a<p)R. Since the quantity 
a<j)Ru (where Ru is the background scalar curvature of the universe) is asymptotically ap- 
proaching A, that means that (1Q7tGn)~^ = l — acf) is asymptotically approaching 1 — A/Rjj, 
which is exponentially large at present (~ 10 120 ). In other words, Newton's "constant" turns 
off in a drastic way. 

This is the problem that seems to afflict any attempt to construct a feedback mechanism 
through a coupling of the cosmon to the scalar curvature through the cosmon potential. 
(This is the conclusion that was reached in the papers and repeated in Weinberg's 

1989 review jjj.) This led to an attempt in 3] to achieve a viable feedback through the 
kinetic term of the cosmon. What is required is that negative powers of R appear in the 
cosmon kinetic term, which, of course, is a peculiar thing indeed jjj. Nevertheless, a model 
that erases the cosmological constant without any fine-tuning of parameters at the classical 
level can easily be constructed as shown in Qj. In the paper we correct and extend the 
analysis in that earlier paper, and address several important issues that were not considered 
there. 

This paper is organized as follows. In section 2 we present the toy model that we are 
discussing and find cosmological solutions for the cases where cosmon energy is dominant 
and where radiation is dominant. (In 3] only the cosmon-dominated case was considered, 
and the analysis of it was faulty, since certain assumptions were made that were not self- 
consistent.) We also discuss the stability of the solutions, and what the cosmon equation 
of state looks like. We find that while there are modifications to the effective Newton's 
constant (because of the appearance of R in the cosmon kinetic term), they are of order 1 
and harmless. 

In section 3, we look at what happens when baryonic matter (or other matter whose 
stress-energy tensor is not traceless) is present. We confirm the conclusion of j^] that this 
does not disturb the feedback mechanism for erasing the cosmological constant. However, 
there is an interesting effect: the cosmon field does wipe out the trace of the stress-energy 
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tensor inside baryonic matter (contrary to what was claimed in [2j]). That is, inside a lump 
of baryonic matter (e.g. a nucleon, atomic nucleus, or neutron star) the cosmon field has a 
pressure and density such that the total stress-energy tensor is traceless. However, this has 
less dramatic effects than one might suppose. This will be analyzed in section 3. 

In section 4, we discuss other issues, such as stability of solutions, what happens if the 
vacuum energy is suddenly reset by a phase transition, and quantum corrections to the 
model and fine-tuning. 

In the Appendix we pent ont certain mistakes and decencies in the ana ly sis in Q. 

2. THE TOY MODEL 

The model that we shall discuss in the rest of this paper has the action 

c 2 

£ = --(<9 CT 0«9»ir 2 + «0 - R - A + Cmatter, (1) 

where c is a parameter with dimensions of M 4 that we take to be of order one in Planckian 
units. Ao is the vacuum energy that we are trying to erase. (It may change suddenly if there 
is a first order phase transition.) The third term is just the usual Einstein-Hilbert action. 
Our metric signature is ( — h ++)• We use the convention that R is negative for a sphere. 

The idea is the following. As the cosmon rolls down its potential hill, the effective vacuum 
energy A = Ao — a<fi will decrease, causing —R to decrease also. This in turn will cause the 
kinetic term of the cosmon to stiffen, so that the rolling of the cosmon will slow down. As 
R approaches zero, so does the velocity of the cosmon. In fact, even though R approaches 
zero, approaches zero faster, so that the kinetic energy of becomes smaller, not larger, 
even though R is in the denominator. Our explicit solutions will show this. 

Having a negative power of R in the cosmon kinetic term [l| can be avoided if we use an 
auxiliary field: 

£ = -lx 2 (d a (j)d a (f)) +a0 

+(d a Xd a X)/X m - bR(XR - c) 2 (2) 

— R — Ao + C-matter- 

The power m in the kinetic term of the auxiliary field X must be large enough so that X 
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can adjust rapidly to changes in R. For large b and m this version of the model reduces to 
that in Eq. (1). We shall therefore use Eq. (1). 
The gravity equation coming from Eq. (1) is 

- /) = -\Tttter + 

- ^(d x (j)d p (t)-\d a (t)d a (t)g xp )R- 2 (3) 
+ (0 A 'D - V( A V>)/, 

where 

/ = c 2 d°<pd a <PR~\ (4) 

and 

A ee A — a<f), (5) 

which implies, of course, that A = —a<p. The last term in Eq. (3) arises from 
<sff~M (~ ^dv4>d CT 'pR ' 2 ^ — '^9Xp)f after two integrations by parts. As we shall see, 
when there is only cosmon energy, the combination of fields that we have called / ends up 
being simply a constant, so that the last term in Eq. (3) vanishes. When there is matter 
present, however, / is not constant and the last term in Eq. (3) becomes very important, 
as we shall see. 

If we take the trace of Eq. (3), we get 

-R(l + §/) = -\T matter + 2A + 3D/. (6) 
The equation of motion for <fi is 

~ --d, [V^(c 2 d^R~ 2 )] = a, (7) 



and assuming a flat Friedmann- Robert son Walker (FRW) ansatz for the metric, the expres- 
sion for the scalar curvature in terms of the scale factor of the universe a(t) is, with our sign 
conventions, 

/- / . \ 2\ 
In I n \ \ 

(8) 




We will solve these equations first in the case where there is no ordinary matter, but 
only cosmon energy. Let us assume that a(t) oc t n and that the fields are constant in space. 
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Then the cosmon equation of motion becomes t 3n d/dt t 3n XR 2 = —a 2 /c 2 , and has the 
solution 

\R-z = -^(a 2 /c?)t-kt-^. (9) 

The second term on the right is a transient that dies off rapidly. (It has already implicitly 
been assumed to be negligible, since its presence would cause the scale factor to deviate from 
a pure power dependence on t.) It is easily shown that the solution to the coupled equations 
gives R proportional to A, so we will write —R = 27A, with the number 7 to be determined. 
Substituting this into Eq. (9) and dropping the transient term gives A = 3 r^r(c 2 /a 2 )t~ 2 , for 
large t Q]. From this, it is easy to see that / is independent oft, so that the last term in Eq. 
(6) may be dropped. Since we are assuming that there is no matter present, that equation 
then takes the simple form —R(l + 1/) = 2A, which implies that (1 + 1/) = 7 -1 . We can get 
another relation between / and 7 from the definition of /, Eq. (4), which we can write as 
/ = (27)~ 3 (c 2 /a 2 )A 2 A~ 3 . Substituting our solution for A into this gives / = (3n + 1)~ 1 7~ 1 . 
Combining the two expressions relating / and 7 gives 7 = 2 ( 6 3n+i) anc ^ / = 6n^r Therefore 

A = lS#Si- 2 . do) 



It remains to use Eq. (8), which tells us that —R = 12n(n — \)t~ 2 . Our solution for A 

2< 6n-i' > o^~ 2- Combining these gives tl 

parameters of the model: 



and 7 gives — R = 27 A = ^g"^ 2 . Combining these gives the power n in terms of the 



6n(6n — l)(n — |) c 2 



(11) 



(3n+l) 2 a 2 ' 

Note that physical quantities must depend on the combination c 2 /a 2 since it is invariant 
under rescalings of the fields. Let us call this combination e. If we take e <C 1, then 



n*± + %e. (12) 



Note that the universe expands faster than t 1 ^ 2 . If there is a small amount of radiation in this 
cosmon-dominated universe, it will redshift as p ra d ~ t~ 2 ~^ e , whereas the "effective vacuum 
energy", A, goes as t~ 2 . We shall find that the same is true in a radiation dominated universe 
with a small amount of cosmon energy: the cosmon/radiation ratio grows as a power of t, 
where the power is of order e. Thus if e ~ 10~ 2 , say, the cosmon energy tracks the radiation 
density very well over many orders of magnitude in time. 

The quantity A contains the vacuum energy A (i.e. the underlying cosmological constant) 
and the cosmon potential energy. If one considers the full "dark" stress energy, which 



includes the contributions from Ao and from both the kinetic and potential energy of the 
cosmon, one finds that for e < 1, p = p/3. In other words, all the "dark" energy - 
cosmological constant plus cosmon — has an equation of state similar to radiation. 

What can we say about Newton's constant? The coefficient of the Einstein tensor in Eq. 
(3), which we may call (16ttGn)~^, is given by (1 — /) = 6(n — \)/ (6n — 1), which for small e 
is approximately ^e. What matters, of course, is not the particular value of (Gjv) e //, which 
is just a meaningless rescaling of the Planck scale, but how (G^eff changes with time. We 
see that if only cosmon energy (including A) is present (Gjv) e // is constant. We shall see 
shortly that when matter is present, (G^eff changes over cosmological time scales, but very 
little if e 1, which is the case of interest. 

Now let us turn to the case of a radiation-dominated universe, with some cosmon energy. 
Since radiation is dominant, one expects the scale factor to grow approximately as t 1 / 2 , so 
we take 

a{t) = a (t/t ) 1/2 (l + gi(t/t ) d + g 2 (t/t ) 2d + ...), (13) 

where to is some time that, as we shall see, is close to when the cosmon energy is equal to 
the radiation energy. The power d is negative, as is gx, so that the cosmon/radiation ratio 
decreases. However, since d turns out to be of 0(e), the cosmon/radiation ratio will change 
very slowly for small e. Eq. (13) gives 

-r = 6d(i + d) 9l t 2 (t/t y 2+d + o(r 2+2d ). (14) 

It is easy to show that the leading behavior of A is t~ 2+2d , so we write 

A = \ t^ 2 (t/t )- 2+2d + 0(t- 2+3d ). (15) 

Substituting Eqs. (14) and (15) into the cosmon equation of motion, Eq. (7), we find 
A = fi ^d~ 9i- Substituting Eqs. (14) and (15) into the Eq. (4) (the definition of /), 
we find / = §§d(l + d)(a 2 /c 2 )gi(t/t ) d to leading order. 

Considering now the trace of the gravity equation, Eq. (6), we see that the leading terms 
on both sides go as t~ 2+d . If we match only the leading terms (which is good enough for our 
purposes), we may neglect A, which contains only terms of order t~ 2+2d and higher. We may 
also neglect the §/ on the left-hand side for the same reason. Since the matter is assumed 
to be radiation, T matter = 0. With fields assumed to be constant in space, we need keep only 
time derivatives, so that Eq. (6) reduces finally to the simple form —R = — 3(J^r + J^)/ = 
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-3d(d + §) ft~ 2 , which gives -R = -§d 2 (d + |)(1 + d)(a 2 / c 2 ) gi t^ 2 (t/t )- 2+d . Combining 
this with Eq. (14) yields 

-±d(l + 2d) = ^ = e. (16) 

For small e, d = — ye, so that 

25 

A/ Pradiation ~ £ 3 € - (17) 

In fact, since the leading term in Einstein's equations gives p ra diation — \t~ 2 , one has 
A/ ' Pradiation = A (t/to) M The preceding calculations give the result A = - ™^% d) gl- 
Consequently, for small e, one has \j Pradiation — ^ e 9i(t/to) 2d - We see, then, that when 
A is comparable to the radiation density, the terms in a(t) of higher order in t d are of the 
same size as the lower order terms, and the whole approximation scheme we have been using 
breaks down, as one would indeed have expected. 

Note that the effective Newton's Constant in an era where radiation dominates is (for 
small e) given by (16ttGn)~L = 1 — f = 1 + 4g 1 (t/t ) _ ^ e . That is, Gn decreases very slowly, 
with Gn/Gn ~ 0(e) ft (although see the discussion in Sec. 4 below). 



3. WHAT HAPPENS INSIDE MATTER 



The mechanism we have been describing relies on the fact that the scalar curvature is 
a "proxy" for the vacuum energy. However, inside a matter distribution there is also a 
contribution from the trace of the matter stress-energy, as shown in Eq. (6). One might 
worry that this would interfere with the mechanism for cancelling the cosmological constant. 

To study this situation, let us consider a spherical mass of radius r m , and constant 
density p m , surrounded by a vacuum. We take p m to be nuclear density, for specificity, and 
the pressure to be negligible. This mass could be a baryon, or a neutron star, for example. 
One might suppose that since in natural gravitational units, p m ~ 10 -80 while A ~ 10~ 120 , 
that p m would dominate the right-hand side of Eq. (6). Let us make this assumption to 
start with, even though it will prove to be false, because it will be instructive and point us 
toward the right solution. 

Using spherical symmetry, we can write the equation of motion for the cosmon inside the 
mass as ^(r 2 §£) = \p 2 m (a 2 /c 2 ). This has the solution, A = A mc + ±p 2 m (a 2 /c 2 )(r 2 - r 2 n ), 
for r < r m , where A mc is the value of A in the vacuum, which we have given in Eq. (10). 
Then AA = A(r = 0) — A mc ~ Pm r m- If r m is the radius of a baryon, this is of order 1CT 120 
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in natural units. Even if r m is the radius of a neutron star, AA is only of order 10~ 80 . Now 
consider the quantity /. Inside the sphere of matter, / = ^(^) 2 R~ 2 ~ p^r^, which for 
a baryon is O(10 -120 ) and for a neutron star is O(10 -80 ). If our assumptions were correct, 
then, it would imply that / = 1 far from the matter and / < 1 inside. But this cannot 
be the case, for then somewhere inside or near the matter distribution the quantity V 2 / 
would have to be of order 1/r 2 ^. For a baryon this is O(10~ 40 ) which is many orders of 
magnitude larger than nuclear density. That is, the last term in Eq. (6) would greatly 
dominate, contrary to our assumption. Thus the assumption that p m dominates is false. 

This shows what the real solution of the coupled equations must look like. It must be that 
/ = fvac + f(x^), where / -C f vac <C 1 and V 2 / = — \p m (inside matter) so that the terms 
— \Tmatter and 3Df on the right-hand side of Eq. (6) cancel. Then —R can have the same 
value inside and outside the matter distribution, allowing A to be nearly spatially constant. 
Consequently, / = f vac even inside matter. This is the only way to get a self-consistent 
solution. 

The net result is the following. The effective Newton's constant (\6tzG^)~^ = 1 — / 
is the same everywhere, both inside baryonic matter (or any other kind of matter) and 
in the vacuum. This is what we would expect anyway from the Equivalence Principle. 
However, something very strange happens; namely the total stress energy tensor acts like 
radiation, since the sum of the matter stress energy and the last term in Eq. (3), namely 
gives something traceless. It is easy to see that if the mass of a ball of 
baryonic matter is M due to T^f atter , the total mass including the last term in Eq. (6) is 
only \M. 

Isn't it catastrophic to have the total stress energy be traceless, even for ordinary matter? 
For example, consider a box of gas. The spatial average of the pressure inside the box would 
be \pgas-, where p gas includes the rest energy of the gas molecules. However, it must be 
realized that this average pressure is almost all coming from pressure in the cosmon field 
inside the material particles. The actual gas pressure due to particles colliding with each 
other is just what would be given by the kinetic theory of gases. On the other hand, the 
box of gas would gravitate as though it had traceless stress energy. 
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4. OTHER ISSUES 



The cosmological solutions we have found are stable under small perturbations. This can 
be seen from the solution to the cosmon equation of motion given in Eqs. (9) and (10). 
The first integral of the equation of motion depends on an integration constant denoted k 
in Eq. (9). One sees that k describes a transient term that dies away with increasing t. 
The integration of Eq. (9) gives another integration constant that we have not shown in Eq. 
(10). However, the exact solution approaches the one given in Eq. (10) for large t. 

There is a simple physical reason why the solutions we have given are stable under small 
perturbations. Imagine that A finds itself at some t with a somewhat smaller value than 
that given in Eq. (10). Then R will also be somewhat smaller than in our solutions. That 
will cause the kinetic term of the cosmon to be more "stiff" , which causes the cosmon field 
to roll down its potential hill more slowly than it does in our solution. In other words, if A 
is initially "too small", it compensates by decreasing more slowly, until eventually it rejoins 
the asymptotic solution given in Eq. (10). In the same way if A starts off somewhat larger 
than the solution in Eq.(10), then R is larger and the cosmon rolls down its hill faster until 
it rejoins the asymptotic solution. 

This feature has an important consequence that would be helpful in making a realistic 
model of cosmology. In a realistic particle theory there would be phase transitions that 
change the vacuum energy density relatively suddenly, such as those connected with the 
breaking of grand unified symmetries, supersymmetry, the weak interaction symmetries, or 
the chiral symmetries of the quarks. Thus, the cosmological "constant" A that has to be 
cancelled off by any cosmon field undergoes sudden jumps downward. The stability of our 
solutions allows our mechanism to accomodate these jumps. Suppose that A = Ao — a<f) 
is evolving as shown in Eq. (10), namely as t~ 2 , and suddenly a phase transition resets 
A to be smaller by an amount of order T 4 , where T is the temperature at which the 
phase transition happens. Now A is of order p ra diation, as we have seen; and p ra diation is of 
order T 4 . So A will be reset downward by an amount that is of order its own size. Let us 
suppose that it is not driven to a negative value by this, but only to a smaller positive value. 
Then, for the reasons we have just explained, the rolling of the cosmon down its hill will 
be slowed until A rejoins the asymptotic solution we have found. In fact, this could happen 
repeatedly at various epochs in the history of the universe: the effective vacuum energy could 
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FIG. 1: Schematic diagram of the evolution of the cosmon field in the presence of phase transitions. 
The solid line represents the actual evolution of the cosmon field while the dashed line is the 
trajectory it would assume in the absence of any phase transitions. In the above case, we have 
included only two phase transitions; after both of which, the vev of the cosmon field asymptotes 
back to the pre-phase transition trajectory. 

be repeatedly reset by phase transitions, and then after a few e-foldings in time, rejoin its 
asymptotic solution. Fig^ schematically illustrates the evolution of just such a cosmon field 
that undergoes several phase transitions. The fact that A levels off for a period of time after 
a phase transition and acts more like a cosmological constant may be a possible explanation 
of the observed dark energy if the universe has recently undergone a phase transition. There 
is a complication that such an explanation would face, namely that the radiation released 
in the phase transition would tend to dominate over A for a period of time. How A and w 
would behave after a phase transition is something that has to be more carefully studied 
numerically before one can tell if this is a viable explanation of dark energy. 

Another effect of the phase transition would be a transient time dependence of Gjy. 
Suppose that A and therefore R is decreased suddenly (compared to t), then by the equation 
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of motion, Eq. one expects <pR~ 2 to be approximately constant during this sudden decrease 
in R. Thus, / = —c 2 (p 2 R~ 3 is driven to a smaller value than the asymptotic solution. As 
A evolves to a rejoin the asymptotic solution, / will be increasing and hence so will Gjq. 
In other words, if we are in an era where A is "recovering" from a phase transition (which 
might explain the near constancy of A), then Gn might be slowly increasing. This can have 
implications for Big Bang Nucleosynthesis. 

It should be noted that if a phase transition were to drive A negative, there would be 
a catastrophe: A could not ever become positive again. Rather it would roll to ever more 
negative values, which would make \R\ larger and the cosmon kinetic term less stiff, so that 
the rolling would accelerate and A would plunge ever faster into the abyss. However, there 
does not seem to be any reason why phase transitions would necessarily have to drive A 
negative. Whether or not one did so would depend on the parameters of the theory. 

There are of course many other issues that would have to be examined to see if a realistic 
cosmology could emerge from the toy model we have described or some similar model. We 
intend to return to these in the future. Another very important issue is whether the kind 
of theory we are describing could make sense at the quantum level. One obvious issue is 
whether there is any fine-tuning is involved in the choice of cosmon kinetic term we have 
made. If one considers the effects of either matter loops or cosmon loops, these do not seem 
to raise any problems of fine tuning (though we have not proved this). One thing that helps 
in this regard is that the cosmon only couples to gravity. However, simple power-counting 
arguments suggest that graviton loops induce a term of the form (d<f)) 2 /R 3 in addition to 
(d<p) 2 /R 2 . The former term would overwhelm the latter unless its coefficient were tuned to 
be less than about 10 -120 , i.e. about the same tuning that is involved in the cosmological 
constant problem itself. However, even if this is the case, that does not mean that the 
model achieves nothing. It seems to be a solution to the cosmological constant problem 
at the classical level (whether a realistic one or not, that remains to be seen). That is, it 
can solve the "small number problem" even if it does not solve the "fine tuning problem" . 
However, any conclusions about the quantum theory are premature, as we have not studied 
it in detail. 
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5. CONCLUSIONS 



We have analyzed a toy model that realizes the idea of a feedback mechanism for cancelling 
the cosmological constant at the classical level. The model has several appealing features. 
First, the total "dark energy" of vacuum energy plus cosmon field energy is always of the 
same order as the energy in radiation and matter. This is certainly suggestive, given the 
observed "cosmic coincidence" that the dark energy is currently of the same order as the 
energy in particles. Second, the mechanism works even if there are several phase transitions 
throughout the history of the universe that reset the cosmological constant. As explained 
in the introduction, any successful cancellation mechanism for the cosmological constant 
must have these features. There are also two interesting modifications of gravity, (a) The 
effective Newton's constant changes slightly over cosmological time scales during periods 
when cosmon energy is subdominant. (b) The cosmon stress energy acts to make the total 
stress energy tensor approximately traceless everywhere, even inside baryonic matter. There 
are several aspects of the model that are less attractive. The most serious is that the scalar 
curvature appears (at least effectively) in the denominator of the cosmon kinetic term. This 
makes it questionable whether this model would make sense at the quantum level. Moreover, 
naive power counting suggests that the model must be fine tuned at the one loop level by 
the same amount that is involved in solving the cosmological constant "by hand" . There is 
also the difficulty that the rolling of the cosmon field tends to erase any dark energy that 
has w = — 1 just as it erases a true cosmological constant. In other words, the cancellation 
mechanism works it too well! However, after phase transitions, the cosmon field tends to stop 
rolling for a while, so that the "effective cosmological constant" is approximately constant 
and may explain the observed dark energy if the universe has recently undergone a phase 
transition. This deserves further study. 

There are several other aspects of the model that also require further study. It must 
be seen what the implications of the tracelessness of the stress-energy tensor are. The 
cosmology of the model must be studied to see how close to being realistic it can be made. 
And whether the model makes sense as a quantum theory needs to be investigated. Finally, 
it would be interesting to see if the feedback idea can be made to work in some other way 
that does not involve such an exotic cosmon kinetic term. 
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APPENDIX 

The present paper is a reconsideration of ideas first presented in a 1987 paper of Barr 
The purpose of this appendix is to point out a number of flaws in the analysis in that 
paper. Before discussing them, one should note that in most of that paper, an auxiliary 
field X was used as in Eq. (2) here. That is, instead of the cosmon kinetic term depending 
on negative powers of the scalar curvature R, it depended on positive powers of X, where 
some potential V(X,R) forced X to be of order 1/R. The cosmon kinetic term assumed 
there was of the form ^(dx^YX 71 , where rj was allowed to be a free parameter for much of 
the analysis, whereas in this paper we always assume that it is 2 (in fact, it can be shown 
that it must be 2 for the power-law solution to dominate). 

(1) The most serious mistake in the analysis in P] was to assume that the scale factor of 
the universe went as a(t) ~ t 1 ^ 2 . In general, as we see from Eqs. (12) and (13) here, this 
is not the case, except in the (unrealizable) limit c 2 /a 2 — > 0. The false assumption that 
a ~ t 1 / 2 together with the assumption that V(X, R) is of the form f(XR), led to the result 
(in Eq. (A15) of Q) that (l&irG N )-} f = |§^y, and (in Eq. (A10) of Q) that the effective 
cosmological constant A ~ t~ 2 '^ r, ~ 1 \ So, in order to have positive Newton's constant, it had 
to be assumed in Q that r] is larger than 2, implying that A falls off more slowly than the 
radiation density (which under the assumption that a ~ t 1 / 2 obviously goes as t~ 2 ). Indeed, 
if rj = 2 + 6, with 6 < 1, then (WttG^j = 0(5) and A/p radiaiion , ~ t°® . 

Now, even though the analysis of |2| was based on a wrong assumption about a(t), and 
the conclusions therefore are invalid, some of those conclusions resemble the results of the 
correct analysis presented here. In the present paper, we also find that A falls off slower 
than the radiation density. And if our parameter e = c 2 /a 2 <C 1, then (IGttGn)^// = 0(e) 

and X/Pradiation ~ t 0( - e \ 
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In P], it was thought that to have positive Newton's constant in the case t] = 2, a rather 
complicated form of V(X,R) had to be used (see section III of that paper). However, that 
is not the case, as we have found here. 

(2) Another mistake in the analysis of Appendix A of Q] was in assuming that V(X, R) 
could have the form f(XR) (cf. Eq. (A2) of [2I). This leads to an inconsistency in the 
analysis as one sees from Eq. (All) of that paper, which is dV/dX = ^g^rjyA 2 . The 
left-hand side is f'(XR)R = O(R), whereas the right-hand side is 0(A 2 ) = 0(R 2 ). That is 
why in Eq. (2) of the present paper we chose the form V(X, R) = bR{XR — c) 2 , which is of 
the form Rf(XR) and does not lead to this inconsistency. 

(3) There is a mistake in Eq. (A3) of Q, in the last term, which contains the operator 
2(g xp \3 — V A V P ). The factor of 2 should not be present, and the indices on V A V P should be 
symmetrized. In other words, the operator should be the same as in Eq. (3) of this paper. 

(4) The analysis of what happens inside matter was incorrect in that it assumed that the 
matter density dominates the right-hand side of what is Eq. (6). As a result, it was thought 
that what we call / here (which corresponds to dV/ dR in |2j if one makes use of the equation 
of motion for X) varies from being very small inside matter to being of order one outside. 
As was noted in j^j] , this would make the terms with second derivatives of / (or equivalently 
dV/dR) enormously large, so large as to dominate over the matter density. As we found 
here, the resolution of this apparent inconsistency is that / inside matter is almost exactly 
the same as outside, but its second derivatives almost exactly cancel the trace of the stress 
energy tensor of the matter. This was not realized in so that some of the discussion in 
section IV of that paper is incorrect. 

In spite of these flaws, most of the major conclusions of that paper stand as correct, 
including the following: (i) in the model where R is contained in the cosmon potential 
Newton's constant is driven to zero; (ii) this problem does not arise if R couples to the 
cosmon through the denominator of its kinetic term; (iii) A falls off slower than the radiation 
or matter density if A dominates; (iv) the feedback mechanism works even when lumps of 
matter are present because the cosmon field does not significantly "sag" inside matter, and is 
therefore everywhere approximately equal to its value in empty space. Our present analysis 
shows that some apparent problems discussed in J^J are not real. Most importantly, the 
disastrously large contributions to the stress energy coming from the second derivatives of / 
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(or dV/dR) do not arise, and there is no difficulty in obtaining a positive Newton's constant. 



[1] For an excellent review, see S. Weinberg, Rev. Mod. Phys. 61, 1 (1989). 
[2] S. M. Barr, Phys. Rev. D 36, 1691 (1987). 
[3] L. H. Ford, Phys. Rev. D 35, 2339 (1987). 

[4] A. D. Dolgov, In "Cambridge 1982, Proceedings, The Very Early Universe", 449-458; T. Banks, 
Nucl. Phys. B 249, 332 (1985); L. F. Abbott, Phys. Lett. B 150, 427 (1985); 

[5] F. Wilczek and A. Zee, unpublished (1983), for more details see the references of [jj; 
F. Wilczek, Phys. Rept. 104, 143 (1984). S. M. Barr and D. Hochberg, Phys. Lett. B 211, 49 
(1988); R. D. Peccei, J. Sola and C. Wetterich, Phys. Lett. B 195, 183 (1987); T. P. Singh 
and T. Padmanabhan, Int. J. Mod. Phys. A 3, 1593 (1988). For more recent works, see 

[6] A. Hebecker and C. Wetterich, Phys. Rev. Lett. 85, 3339 (2000) |arXiv:hep-ph/0003287| ; 
M. Sami and T. Padmanabhan, Phys. Rev. D 67, 083509 (2003) [Erratum-ibid. D 67, 109901 
(2003)] |arXiv:hep-th/0212317| ; S. Nojiri and S. D. Odintsov, Phys. Lett. B 599, 137 (2004) 
|arXiv:astro^ph7 0403622 ; Phys. Lett. B 599, 137 (2004) |arXiv:astro-ph/0403622"1 ; I. Dym- 
nikova and M. Khlopov, Grav. Cosmol. Suppl. 4, 50 (1998); Eur. Phys. J. C 20, 139 (2001); 
Mod. Phys. Lett. A 15, 2305 (2000) |arXiv:astro-ph/0102094| .; M. Pavsic, Found. Phys. 35, 



1617 (2005) arXiv:hep-th/0501222 ; Phys. Lett. A 254, 119 (1999) |arXiv:hep- th/9812123 . 



[7] S. R. Coleman and F. De Luccia, Phys. Rev. D 21, 3305 (1980). 



[8] T. Banks, arXiv:hep-th/0211160 



[9] For papers where negative powers of R appear in the Lagrangian, see S. M. Carroll, V. Duvvuri, 
M. Trodden and M. S. Turner, Phys. Rev. D 70, 043528 (2004) |arXiv:astro-ph/0306438"] ; 
D. N. Vollick, Phys. Rev. D 68, 063510 (2003) |arXiv:astro-ph/0306630| ; T. Chiba, Phys. 
Lett. B 575, 1 (2003) |arXiv:astro-ph/0307338| ; J. M. Cline, |arXiv:astro-ph/0307423 



S. Nojiri and S. D. Odintsov, Unification of the inflation and of the cosmic accelera- 
tion," Phys. Rev. D 68, 123512 (2003) |arXiv:hep-th/0307288|; X. Meng and P. Wang, 



Class. Quant. Grav. 20, 4949 (2003) |arXiv:astro -ph/0307354| ; Class. Quant. Grav. 21, 951 



(2004) |arXiv:astro-ph/0 308031 ; A. D. Dolgov and M. Kawasaki, Phys. Lett. B 573, 1 
(2003) |arXiv: astro-ph/0307285| and E. E. Flanagan, Phys. Rev. Lett. 92, 071101 (2004) 
|arXiv:astro-ph/0308m"] . 



17 



[10] Similar time dependence in another kind of model was found in O. Bertolami, Nuovo Cim. B 
93, 36 (1986). 



18 



